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We consider a higher order term in the 5N expansion for the CMB power asymmetry generated 
by a superhorizon isocurvature field fluctuation. The term can generate the asymmetry without 
requiring a large value of /nl- Instead it produces a non-zero value of gNL- A combination of 
constraints leads to an allowed region in /nl — (?nl space. To produce the asymmetry with this 
term without a large value of /nl we find that the isocurvature field needs to contribute less than 
the inflaton towards the power spectrum of the curvature perturbation. 


I. INTRODUCTION 

Inflation is widely accepted as the likely origin for 
structure in our universe. Its generic predictions of a 
nearly scale invariant and close to Gaussian primordial 
curvature perturbation, C,, have been confirmed with in¬ 
creasing precision by successive Cosmic Microwave Back¬ 
ground (CMB) experiments. There are, however, also ob¬ 
servational anomalies which are harder to explain within 
the standard inflationary paradigm. One such anomaly 
is the hemispherical power asymmetry - the observa¬ 
tion that for scales with I < 60 there is more power in 
CMB temperature fluctuations in one half of the sky than 
the other. First identified in the Wilkinson Microwave 
Anisotropy Probe data [DU, it was later confirmed by 
the Planck collaboration [5] and others |6H8], although its 
significance remains disputed [^. In this work we treat 
the asymmetry as a real effect which requires a primor¬ 
dial origin. So far, CMB data has been fitted to a tem¬ 
plate which models the asymmetry as a spatially linear 
modulation. 

The leading primordial explanation for this asymmetry 
is the Erickcek-Kamionkowski-Carrol (EKC) mechanism 
[Mill], in which a long-wavelength isocurvature pertur¬ 
bation modulates the power on shorter scales. Further 
work investigating this effect includes Refs. [ISHM- The 
origin of the long wavelength mode may be explicitly re¬ 
alised in the open inflation scenario of m or due to a 
domain wall, as in, for example, I1I1I21. 

The 6N formalism provides a convenient expression 
for the modulation of power by a super-horizon mode, as 
reviewed below. In principle many terms in this SN ex¬ 
pansion can contribute to the observed asymmetry. Until 
now, however, most theoretical work has focused on the 
leading term, which can have the form of a spatially lin¬ 
ear modulation. 

If the leading term in SN is responsible for the asym¬ 
metry then a further consequence is that the local bis¬ 
pectrum parameter must satisfy the constraint /nl ^ 30/? 
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m m on the scales that are modulated. A value of 
/? < 1 can be achieved but only if our observable uni¬ 
verse is located at a fine-tuned region within the long- 
wavelength perturbation |12] . and otherwise can be much 
larger than one. Combined temperature and polarization 
data bounds a purely scale-independent local bispectrum 
as /nl = 0.8 ± 5.0 at 68%CL |24], while we work with 
/nl < 10 as a rough 95%CL. The asymmetry appears to 
be scale dependent jSS], and hence the non-Gaussianity 
produced must also be, but there are no direct constraints 
on such a strongly scale dependent non-Gaussianity. A 
new parametrisation of the scale-dependence of the non- 
Gaussianity and its application to the scale-dependence 
of the asymmetry was given in |26| . which includes an 
accompanying ^nl- It is, however, perhaps unlikely that 
a very large value of /nl could be accommodated by cur¬ 
rent observations, even if /nl decays with scale. 

In this short paper, therefore, we investigate whether 
the next term in the SN expression for the asymmetry 
could instead be responsible. We And it can, without 
violating any other observational or self-consistency con¬ 
straints. It contributes a more general modulation of the 
power, leading to an asymmetry, which does not necessar¬ 
ily only involve a spatially linear modulation m- Using 
this higher order term requires a non-zero value of 5 nl, 
but allows for a smaller value of /n l than when the linear 
term alone contributes. If this higher order term is re¬ 
sponsible for the asymmetry, then the allowed parameter 
space indicates the modulating isocurvature field must 
contribute less than the inflaton towards the total power 
spectrum of the curvature perturbation on scales which 
are modulated, and this may be considered a fine-tuning 
of the model. Related to this, we And that if this higher 
order term is dominant in our observable patch, then in 
certain neighbouring patches the linear term will instead 
be dominant. 

In this paper, as a first step we only focus on one of 
the higher order terms, but the idea is more general and 
could be applied to a combination of higher order terms. 
Satisfying the constraints in that case might be more 
complicated than the simple use of exclusion plots that 
we employ here. 
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II. GENERATING THE ASYMMETRY 

A. The SN Formalism 

Our calculation is performed within the 6N formalism 
[2SH35] which states that ^ can be associated with the 
difference in the number of e-folds undergone by neigh¬ 
bouring positions in the universe from an initial flat hy¬ 
persurface at horizon crossing to a final uniform density 
one when the dynamics have become adiabatic: C = 

On the flat hypersurface the inflationary fields are not 
constant, and by writing as a function of the fields, 
6N can be written as a Taylor expansion in the horizon 
crossing field fluctuations. 

We consider two scalar fields, though our work easily 
generalises for more than two fields, and we take both 
our fields to have canonical kinetic terms. We choose the 
inflaton field, denoted (j), to be the direction in held space 
aligned with the inflationary trajectory at horizon exit, 
t*, so that e* = and this implies the derivative of N 
with respect to the inflaton is a constant 

= (2e*)-i/2. (1) 

The isocurvature held orthogonal to (f) is denoted y, and 
the curvature perturbation has contributions from both 
fields 

C = + ■■■ (2) 

where we have neglected terms with higher order cj) 
derivatives since they are negligible. The arguments of 
N and its derivatives are usually taken to be the aver¬ 
age values of the fields within our observable universe, 
denoted and xoj while S(j) and Sx contain all fluctu¬ 
ations in (/) and x with wavelengths of order the size of 
our observable universe or less. 

The power spectrum of the curvature perturbation is 
then given by 


P(- = NjNj 



( 3 ) 


where / runs over {</>, x}) the summation convention has 
been used, and we have neglected higher order d4> and 6x 
correlators. 

Non-Gaussianities in ( are generated because of the 
non-linear relationship between (j and (5x in ([^- In par¬ 
ticular, one finds for the local bispectrum, /atl, and 
trispectrum, gNL, parameters that [551151] 


5 N^xN^^ 

6 {NjNj^ 


_ 25 N,xxx^^x 
~ 54 {NjNj f ■ 


( 4 ) 

( 5 ) 


In what follows we will only be concerned with the mag¬ 
nitude of fNL and gNL , I/tvlI and but to avoid 


clutter we will drop the absolute symbols. We will also 
use the expression for the tensor-to-scalar ratio 
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NjNj 


( 6 ) 


and we will And it convenient to define the contribution 
of X to the total power spectrum 


X = — 

Pc 


Kx 1 ^ 

NjNj 16e* ■ 


( 7 ) 


B. Superhorizon Fluctuation 

In addition to the background value of the fields inside 
our observable universe, {^OiXo} and their fluctuations 
with wavelength inside our observable universe, {Scjj, Sx}-, 
the EKC mechanism works by postulating a superhorizon 
field fluctuation in x, denoted Ax(x), with wavelength, 
much larger than the size of our observable universe, 
this size given by the distance to the last scattering sur¬ 
face, Xd, such that fc^Xd <C 1. We assume the leading 
order behaviour Ax(x) = Ax(n ■ kL) for x within our 
observable universe, where n = x/|x| and Ul = kL/|kL|, 
and we don’t assume any particular form for the fluctu¬ 
ation outside of our observable patch. Note that in this 
paper we take Ax to be the maximum variation in x 
across our patch about our observable universe’s average 
field value xo as seen in the left panel of Fig [55], 
Superhorizon fluctuations source multipole moments 
in the CMB, upon which there are constraints from the 
observed homogeneity of the universe HOI E]. Using 

the non-linear results of m, together with the multipole 
constraints from m, we have the following homogene¬ 
ity constraints from the quadrupole and octupole respec¬ 
tively |5S] 

(A^)" I < 1.1 X 10-4 (8) 

I <8.6x10-4 (9) 

where we have assumed no cancellation between SN 
terms. We also take the following constraint 

< aPl'^ (10) 

where Pq = 2.2 x 10’ -9 m and a is some threshold pa¬ 

rameter. 


C. Asymmetry 

The superhorizon fluctuation modulates the power 
spectrum on shorter scales, and so it depends on the di¬ 
rection n through 

Pci^] =7’c[Xo + Ax(n)]. 


( 11 ) 
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Since Ax(n) < in our patch, and Ax is small, we can 
Taylor expand in (11) in powers of Ax(n) giving 


Pc[fi] =Pc l + 


m—1 


where the round brackets indicate multiplication. 


= 


1 (Ax)^ 


2P 




9X" 


( 12 ) 


(13) 


and we have used the shorthand that when P(^ and its 
derivatives appear without an argument they are taken to 
be evaluated at the average field values of the observable 
universe. Observations indicate a power asymmetry, with 
the power along the preferred direction ri = kL being 
greater than the power on the opposite side of the sky 
n = —kL- We note that only the odd m terms in (12 1 can 


contribute towards an asymmetry of this sort, with the 
even terms contributing only towards general anisotropy. 

Usually only the m = 1 term is kept, and the data has 
been fitted to this with the result that ^ Ai — 0.07. 
The m = 1 and m = 2 terms were considered in |38| . 
Here we consider instead the m = 3 term, since this can 
contribute towards asymmetry [39]. Ideally a fit to the 
data with m = 1,2,3 terms should be done to constrain 
the parameters Ai, A 2 and A 3 . In the absence of this, we 
will look at the simplest case involving only the m = 3 
term and take |3D] A 3 > 0.07. 


Differentiating ^ three times gives 

^‘',xxx _ ,xxx^,xx + 2A^,xxxxA^.x 
Pc NjNj ■ ^ ’ 

We will be interested in the case where the asymme¬ 
try is generated by the N,xxx^,xx t^rm, and we neglect 
N^xxxx^ so that our asymmetry is given by [35] 


_ ^,xxx^,xxi^x)^ 
^ 2NjNj 


(17) 


In this case, we now show there is still a lower bound on 
Jnl, but this time it depends on x defined in ([^. Using 
together with the octupole constraint (|^, we find 


/nl 


6 {NjNj)^ 


>9.5 



(18) 


We see that if x is sufficiently small, we can have an ac¬ 
ceptably small fixL in this scenario. We will later show 
that there is a lower bound x > A 3 , and so 9.5 is the small¬ 
est value of /ml allowed from this cubic term alone |43| . 
which is an improvement compared to the contribution 
from the linear term alone. 


F. Consistency Checks 


D. Linear Term Asymmetry 


For simplicity we assumed that the asymmetry is only 
due to the m = 3 term in (12 1 , which then must be larger 
than the to = 1 term. We therefore require 


It has been noted in e.g miiiiiTni that a large f^L 
accompanies the asymmetry when only the to = 1 term 
is considered, and we briefly review this now. Differenti¬ 
ating ^ gives 


^ ^,xx^,x^X 
^ NjNj 


We now combine this with constraint ^ giving 


Jnl 


^N.xxKx 

0{NjNjY 


>37 



(14) 


(15) 


which is outside of the observational bounds for a local- 
type non-Gaussianity |41| . 


E. Cubic Term Asymmetry 


The asymmetry may be due to multiple odd to terms 
in (12). We will now show that postulating the cubic 


TO = 3 term is dominant over the linear to = 1 term, and 
is responsible for the asymmetry, allows the constraint 
on fixL to be relaxed, but introduces new ones on gNL- 
Later we will check the self-consistency of ignoring the 
TO = 1 term compared to the to = 3 one. 


-^,xxxA^.xx(Ax)^ 

‘^NxxNx 


> 1 . 


(19) 


Even powers of Ay don’t contribute towards the asym¬ 
metry but they do still cause more general anisotropy of 
the power spectrum. Since these anisotropies have not 
been observed, we also demand the following 


A^.xxx-^.xxAx 

N N 
^',xxx^^ ,x 


> b 


and ^-xxxjV.xxAx ^ ^ 

AXX 


( 20 ) 

( 21 ) 


where b, c are some threshold parameters. 

There is a lower limit on a: = x{xo) coming from x{xo — 
Ay) > 0, by definition ([^. Expanding out a;(xo — Ay) 
to cubic order and neglecting the linear term, we find 
a;(xo)>^3 for b,c^O{l). 


G. Allowed Parameters 


We have six constraints to simultaneously satisfy: (|^, 
([§, (1^, (l2§ and (ig. Using to substitute 
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Figure 1. Left: The value of x varies by an amount Ax 
from its average value, xo, within our observable universe (the 
interior of the two vertical lines), due to the long-wavelength 
fluctuation (solid wave). The average within our observable 
universe (long-dashed line) is not necessarily the same as the 
background value over the entire universe (dotted horizontal 
line). Right: Exclusion plot for (22l-dark-blue, (23horange, 
( |24[ |-green, ( |25[ )-red, (|26[)-purple, ( |27p -brown, with x = 0.07, 
a = 10, b = 0.25 and c = 1. The red and brown lines are hard 
to see on this scale at the bottom of the plot. The allowed 
region is left white. 


for Ax, and using (0, (0) and (§ the six constraints 
become, respectively. 


9nl > 

fNL > 

9nl > 

9nl > 
9nl < 

9nl > 





( 22 ) 

(23) 

(24) 

(25) 

(26) 
(27) 


In the right panel of Fig l] we plot the allowed region, 
left in white, for (22l-(27l, with x = = 0.07, a = 

10, b = 0.25 and c = 1. We find that the cubic term 
can generate the required asymmetry with a lower value 
of fwL than from the linear term alone. Moreover it 
requires a non-zero value of > 5 x 10^ for the smallest 
allowed values of /nl- Note that if x is much bigger than 
0.07 then this pushes the allowed values of /nl and ^nl 
up. The small value of x = 0.07 may be considered a fine- 
tuning required when only the m = 3 term generates the 
asymmetry. 


H. Outside Our Observable Patch 


In the above we neglected the first order m = 1 term in 


( 121 , assuming that this term is small in our observable 
universe. However, since we are considering a scenario in 
which W and are non-zero, neighbouring regions 

of the universe with a different background field value 
may have a larger first order term. This is closely related 
to a similar effect in inhomogeneous non-Gaussianity m- 
I48| . If this term is larger in neighbouring patches this 
would not violate observational bounds, but would imply 
that our position within neighbouring regions was finely 
tuned - in the sense that neighbouring regions would 
instead see a dominant first order term. Although not 
invalidating the proposed scenario, it would make it less 
appealing. The biggest change in the average value of 
X is in a neighbouring patch along the direction of the 
long wavelength mode, where its average value is of order 
Xo + Ax, since Ax > Sx- The first order term in these 
patches is then of order 


^A,A- . A ^AA lA 

Xo+Ax 


.Ax(Xvx + ^.xxxXx) 


+ 2 Ax) XxxxXxx 


Xo 


(28) 


where we have neglected fourth and higher derivatives of 
N. The order Ax term in (28 1 is related to the zeroth 
order term by 


^X(-N,xx XxxxXx) 


XxxXx 


Xo 


> (6 + c) 


(29) 


and so these terms are of comparable order for 6 , c = 0 {1) 
and if (20 1 and (21) are not hierarchical inequalities. The 
order (Ax)^ term in (281 is related to the zeroth order 
term using (191 


3(Ax) XxxxXxx 


^XxxXx 


> 3 


Xo 


(30) 


so we see that the first order term in Ax in ( 12 ) in these 


neighbouring patches will actually be of the same order 
or larger than the cubic one in our own patch which we 
consider to be repsonsible for the asymmetry. This then 
implies that in these neighbouring patches the value of 
fNL is necessarily larger than in our own patch. This 
agrees with the result of j49] that if qnl 'A> Inl Ai our 
observable patch, then neighbouring patches will generi- 
cally have a larger value of /nl than in our own. If the 
asymmetry in our patch is due to the third order term 
rather than the linear term, then our patch should be 
considered fine-tuned compared to its neighbours along 
the direction of the long wavelength mode. 
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III. CONCLUSION 


We have presented a mechanism involving a mod¬ 
ulating isocurvature field which can produce the re¬ 
quired hemispherical power asymmetry while satisfying 
the homogeneity constraints, and which produces non- 
Gaussianity within observational bounds. A novel fea¬ 
ture is the non-zero value of required to generate 
this asymmetry. We note that there are models with a 
large gNL and small Jnl, for example, m and m- A 
requirement on the model is that the isocurvature field 
contributes a small amount towards the power spectrum 
of the curvature perturbation, which could be consid¬ 
ered a fine tuning. We also note that the large minimal 
value of gNL required implies our observable patch of 
the universe has a significantly smaller value of /nl than 
our neighbours. The observed asymmetry is scale depen¬ 
dent, with a smaller asymmetry on small scales, which 
this model does not account for. 

If the observed asymmetry is due to the higher order 
term considered in this work, then this will put strong 
bounds on /nl and gNL- Measurements of /nl and gNL 
outside of our allowed region would falsify models which 
use this cubic term to generate the asymmetry. 

The cubic term has a different fi-dependence com¬ 


pared to the first order term. For this paper we assumed 
A 3 > 0.07, but we would like to see a fit to the data with 
the m = 1,2,3 terms, in order to properly constrain the 
parameters Ai, A 2 and A 3 . 

This study has shown that a higher order term can gen¬ 
erate the required asymmetry, relaxing the constraint on 
/nl compared to that generated only by the first order. 
Perhaps the other cubic order term in (16 1 , A^.xxxx-A.xi 
may also contribute - although the bound on the non¬ 
linear parameter, Hnl |52], associated to this term is 
considerably weaker than that on gNL-, and so this term 
is not as easily falsifiable. Indeed, since the third order 
term can have a large contribution, other higher order 
terms (and combinations of them) may also be signifi¬ 
cant. Our work prompts investigation of the case where 
6 N can’t be Taylor expanded. 
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